The magnetic state of Yb in Kondo-lattice YbNi 2 B 2 C 
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We report neutron scattering experiments performed to investigate the dynamic magnetic prop- 
erties of the Kondo-lattice compound YbNi2B2C. The spectrum of magnetic excitations is found to 
be broad, extending up to at least 150 meV, and contains inelastic peaks centred near 18meV and 
43meV. At low energies we observe quasielastic scattering with a width F = 2.1 meV. The results 
suggest a Yb 3+ ground state with predominantly localized 4/ electrons subject to (i) a crystalline 
electric field (CEF) potential, and (ii) a Kondo interaction, which at low temperatures is about an 
order of magnitude smaller than the CEF interaction. From an analysis of the dynamic magnetic 
response we conclude that the crystalline electric field acting on the Yb ions has a similar anisotropy 
to that in other i?Ni2B2C compounds, but is uniformly enhanced by almost a factor of 2. The static 
and dynamic magnetic properties of YbNi2B2C are found to be reconciled quite well by means of an 
approximation scheme to the Anderson impurity model, and this procedure also indicates that the 
effective Kondo interaction varies with temperature due to the crystal field splitting. We discuss 
the nature of the correlated-electron ground state of YbNi2B2C based on these and other experi- 
mental results, and suggest that this compound might be close to a quantum critical point on the 
non-magnetic side. 



I. INTRODUCTION 

One of the delights of /-electron metals is their ten- 
dency to show beautiful electronic ordering phenomena 
at low temperatures. Typically one can expect to find 
long-range magnetic ordering, superconductivity, heavy- 
electron behaviour, or quantum critical effects, and some- 
times two or more of these in the same material. The 
character of the ground state is often governed by several 
different microscopic interactions of comparable strength. 

The diversity of possible electronic phases is exempli- 
fied by the i?Ni2B2C family, R being Sc, Y, most lan- 
thanides, and some actinidesi2iiii Most members of this 
family exhibit either magnetic ordering or superconduc- 
tivity, and for R = Dy, Ho, Er and Tm there exists 
a temperature range over which superconductivity and 
magnetic order coexist. These four members of the fam- 
ily have been the subject of extensive investigations into 
the interplay between magnetism and superconductivity. 
The materials with R = Sc, Y, Lu, Th and possibly Ce are 
non-magnetic superconductors^^ whereas LaNi2B2C is 



a conventional metal. Among the remaining members of 
the series, only YbNi2B2C displays neither magnetic or- 
dering nor superconductivity. This compound is metallic 
with enhanced electronic transport and thermodynamic 
coefficients at low temperatures, and for this reason it 
has been classified a heavy fermion systemic 

This paper addresses the nature of the ground state of 
YbNi2B2C. Since it is stoichiometric, YbNi2B2C is one of 
a rather small number of known Yb Kondo-lattice com- 
pounds. Its properties are then expected to be strongly 
influenced by local Kondo singlet formation, with com- 
petition from long-range magnetic ordering of localized 
4/ moments coupled by the RKKY exchange interac- 
tion. Given the right energy balance, Kondo-lattice com- 
pounds can exhibit a coherent Fermi liquid phase at low 
temperatures with strongly renormalized quasiparticles, 
or 'heavy fermions'. In the heavy fermion state an abun- 
dance of low-energy spin fluctuations mediate the inter- 
actions between the quasiparticles and strongly influence 
the physical properties. An additional factor is the crys- 
talline electric field (CEF), which by splitting the 4/ lev- 
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els dictates the symmetry of the 4/ ground state and 
makes the effective 4/ magnetic moment temperature de- 
pendent. 

Evidence for a coherent heavy fermion state 
in YDM2B2C has been found in several different 
propertiesA^ The electronic heat capacity has a peak 
near 8 K, and decreases below 5 K with a large 
slope corresponding to a Sommerfeld coefficient 7 « 
0.5 J mol -1 K -2 £ This value of 7 translates into a single- 
impurity Kondo temperature Tk ~ 10 K. The resistiv- 
ity decreases steadily below room temperature and then 
more sharply below a temperature in the range 10-40 K 
depending on the annealing treatment of the sample^. At 
temperatures below ~1.5K the temperature dependence 
of the resistivity exhibits a positive curvature reminiscent 
of the T 2 variation characteristic of a strongly-correlated 
Fermi liquid. The temperature dependence of the spin- 
lattice relaxation rate measured by n B nuclear magnetic 
resonance^ shows a transition from local moment relax- 
ation at high temperatures to Korringa-likc behaviour 
below 5 K suggestive of a Fermi liquid with a high den- 
sity of states. Finally, no magnetic ordering has been 
observed in YbN^I^C down to 0.023 Kjii a tempera- 
ture much less than the Yb magnetic ordering temper- 
ature of 0.4 K predicted by de Gennes scaling from the 
heavy R ions in i?Ni2B2C, and nearly four orders be- 



low the Weiss temperature V 
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Curie- Weiss law fit to the high-temperature susceptibil- 
ity. These observations suggest the existence of a strong 
antifcrromagnetic interaction between Yb and conduc- 
tion band states, which causes a screening of the Yb 4/ 
moment at low temperatures by the Kondo effect. 

Here we describe neutron inelastic scattering measure- 
ments of the magnetic excitation spectrum of YbNi2B2C 
in the 0-200 meV energy range. This technique directly 
measures the spin fluctuation spectrum, and allows us 
to determine the Kondo, CEF and RKKY energy scales 
quantitatively. Our main findings are (i) that the CEF 
potential is larger than either the Kondo or RKKY in- 
teractions, and (ii) that the magnetic excitation spectra 
and other bulk magnetic properties can be quite well un- 
derstood in terms of a model containing CEF and Kondo 
interactions provided that the latter is allowed to increase 
with temperature. A preliminary account of our experi- 
mental data was reported in Ref. 0. Neutron spectro- 
scopic data have also been published by Sierks it et al.fi 3 - 
and though our conclusions are partly in accord with 
those of Sierks et al. there are also differences which we 
will seek to explain later. 



II. SAMPLE PREPARATION 

Experiments were performed on single-crystalline and 
polycrystallinc YbNi2B2C. All the samples were pre- 
pared with boron enriched to 99.5% n B to reduce neu- 
tron absorption by the 10 B isotope present in natural 
boron. Single-crystalline LUM2B2C and polycrystalline 



Y0.5LU0.5N12B2C were used as non-magnetic reference 
samples. 

The single crystals were grown at Ames Laboratory by 
the high-temperature M2B flux methodi^i^ In the case 
of YbNi2B2C, excess Yb was added to compensate for the 
loss due to evaporation. The crystals were plate-like with 
typical dimensions 4 x 4 x 0.3 mm 3 , and the mass of the 
largest crystal was 0.25 g. For the neutron experiments 
we prepared a mosaic of 40 crystals with a total mass of 
approximately I g. The crystals were co-aligned to within 
5° and glued onto a thin sheet of aluminium. A similar 
mosaic of LUM2B2C crystals was used to estimate the 
non-magnetic signal. 

Polycrystalline YbNi2B2C was prepared by a solid- 
state reaction method similar to that described by Dhar 
et al£ and polycrystalline Y0.5LU0.5N12B2C was prepared 
by the standard arc-melting technique under flowing ar- 
gon. X-ray powder diffraction was used to check for im- 
purities, and the YbNi2B2C sample was found to contain 
a measurable amount of Yb2C>3. Subsequently, a quan- 
titative analysis of the composition of this sample was 
carried out by multi-phase refinement of neutron powder 
diffraction data collected on the Dlb diffractometer at 
the Institut Laue-Langevin (ILL). The amount of Yb203 
impurity was found to be 9% by mass. As described be- 
low, we corrected for the spurious signal produced by this 
impurity by performing measurements on commercially- 
obtained Yb2C-3 powder (99.9% purity) under identical 
conditions as used for YbNi2B2C, with Y2O3 powder as 
the non-magnetic reference. 



III. MAGNETIC EXCITATION SPECTRUM 

A. Neutron inelastic scattering cross-section 

The quantity measured by neutron inelastic scattering 
is the differential cross-section per unit solid angle £1 and 
scattered neutron energy Ef, given by 



d 2 a _ kt 
dndEf ~ hi 



(f) 



Here, k\ and kf are the incident and scattered neutron 
wavevectors, S"(Q, u>) is the response function of the sam- 
ple, Q = k ; — kf is the scattering vector, and hid = E{—Ef 
is the energy transferred from the neutron to the sample. 
From linear response theory it can be shown that the dy- 
namical part of the response function for a paramagnetic 
ion is given in the dipole approximation by«i£ 



(7^0) 



cxp{-2IT(Q)}|/(Q)| 2 lu{1 + «(«)} 



- Q 2 a )x aa F aa (uj), 



(2) 



where 7 = — f .913, r = 2.818 x f0~ 15 m is the classical 
electron radius, exp{— 2VT(Q)} is the Debye- Waller fac- 
tor, |/(Q)| 2 is the squared modulus of the magnetic form 
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factor, n(uj) = 1 / {exp(fi,w / k^T) — 1} is the Planck dis- 
tribution, Q a is the a (= x, y, z) component of the unit 
scattering vector, \ aa is a leading-diagonal element of 
the static single- ion susceptibility tensor—, and F aa (uj) 
is a spectral-weight function with unit normalization: 

/oo 
F aa {uj)dLU = 1. (3) 

In the case of a polycrystalline sample, Eq. J5J must be 
averaged over all orientations. If, in addition, we extrap- 
olate the scattering to zero Q then the expression for the 
dynamic part of the response function becomes 

S(fl,w) = ^■«{l + n(a;)}2x.„F(i 1 ;) ) (4) 

where Xav = \ {x xx + X yv + X zz ) ls * ne powder-averaged 
static susceptibility. 

B. Experimental details 

We employed four different neutron spectrometers to 
study the magnetic excitation spectrum. The polycrys- 
talline samples were measured on the high energy transfer 
(HET) chopper spectrometer at the ISIS spallation neu- 
tron source and on the IN5 chopper spectrometer at the 
Institut Laue-Langevin (ILL). These experiments probed 
the energy ranges 5-200 meV and 0.3-2.5 meV respec- 
tively. The single-crystal samples were measured on the 
IN 14 triple-axis spectrometer at the ILL and on the IN6 
time-of-flight spectrometer, also at the ILL. These latter 
experiments provided information on the Q dependence 
of the low-energy excitations in the energy ranges 0.2— 
6meV (IN14) and 0.2-4.5 meV (IN6). All four spectrom- 
eters were equipped with a variable-temperature liquid 
helium cryostat allowing measurements to be made as a 
function of temperature. Experimental details specific to 
each spectrometer are as follows: 

(1) HET. Data were collected in three runs, with neu- 
trons of incident energy 35meV, 75meV and 250 meV 
respectively. Spectra recorded in banks of detectors dis- 
tributed around the incident beam direction were av- 
eraged. The mass of polycrystalline YbNi2B2C in the 
beam was approximately 13 g. To provide an estimate 
of the non-magnetic background scattering we also mea- 
sured spectra at the same incident energies from a similar 
mass of polycrystalline Yo.5Luo.5Ni 2 B 2 C. This composi- 
tion was chosen because it has the same cross-section for 
neutron absorption as YbNi2B2C. 

(2) IN5. Measurements were made with the incident 
energy fixed at 3.1 meV. The data did not show any de- 
tectable Q dependence apart from a very slow reduction 
in intensity with |Q| for a given energy consistent with 
the variation of the magnetic form factor, and so the 
counts recorded in the whole detector bank (extending 
from ^15° to ^130° in scattering angle) were averaged 
to improve statistics. 



(3) IN14. Spectra were recorded by scanning the in- 
cident neutron energy with a fixed final energy of Ef — 
4.7 meV. The incident and final energies were selected by 
Bragg reflection from arrays of pyrolytic graphite crys- 
tals. A beryllium filter was placed immediately after the 
sample to suppress higher-order harmonics in the scat- 
tered beam. Two settings of the crystals were used, giv- 
ing access to the (h, h, I) and (h, 0, 1) planes in reciprocal 
space. 

(4) IN6. The incident neutron energy was fixed at ei- 
ther 3.1 meV or 4.9 meV. As on IN5, we averaged the 
neutron counts recorded in the whole detector (~10° to 
~115°). This means that the recorded spectra corre- 
spond to an average over a range of Q. We chose two 
orientations of the crystal relative to the incident beam 
so that in one case the average Q was approximately par- 
allel to the a axis of the crystal, and in the other it was 
approximately parallel to the c axis. 

On all the time-of-flight spectrometers (HET, IN5 and 
IN6) the scattering from a standard vanadium sample 
was used to normalize the data in different detector banks 
and to convert the spectra into units of absolute scatter- 
ing cross-section. 



C. Results 

We begin with the high-energy data collected on HET. 
Figurenshows an example of the raw data collected with 
an incident energy of 250 meV and a sample temperature 
of 10 K. The data are from a detector bank covering a 
range of scattering angle <fi from 3° to 7°, i.e. ((f)) = 5°. 
Spectra from both YbNi2B2C and Yo.5Luo.5Ni2B2C are 
shown. 

To a first approximation, the difference between the 
scattering from the two samples is the magnetic scatter- 
ing, but a straight subtraction is not an entirely satis- 
factory way to isolate the magnetic scattering because 
(a) the nuclear scattering amplitudes of Yb, Y and Lu 
are different, and (b) there are small differences in the 
phonon spectra of the two materials, particularly at low 
energy. For a more accurate estimate of the non-magnetic 
scattering from YbNi2B2C we adopted the following pro- 
cedure. First we took the ratio of the energy spectra from 
YbNi 2 B 2 C and Y .5Luo.5Ni2B 2 C measured in a high- 
angle (~130°) detector bank. At high angles the spec- 
trum measures the phonon density of states weighted by 
the scattering power of the elements. The magnetic scat- 
tering is negligible because of the decay of the magnetic 
form factor with |Q|. Multiplication of this high-angle ra- 
tio by the low-angle spectrum of Yo.5Luo.sNi2B2C then 
gives a good estimate of the non-magnetic low-angle 
background of YbNi2B2C. This procedure is based on 
the assumption that the non-magnetic inelastic scatter- 
ing at low angles is dominated by (elastic + 1-phonon) 
multiple scattering and therefore closely resembles the 
scattering measured at high angles. The validity of this 
assumption for high-energy neutrons incident normally 
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FIG. 1: Raw neutron scattering energy spectra obtained from 
polycrystalline samples of YbNi 2 B 2 C and Y0.5LU0.5N12B2C 
on the HET time-of-fiight spectrometer. The incident neu- 
tron energy was 250 me V, and the data shown were recorded 
in the HET low-angle detector bank and normalized per for- 
mula unit (f.u.) of YbNi 2 B 2 C. The YbNi 2 B 2 C non-magnetic 
background scattering (filled circles) was estimated by the 
method described in the text. 



FIG. 2: Magnetic response function of polycrystalline Yb 2 03 
measured on the HET time-of-flight spectrometer. The data 
are corrected for (i) the non-magnetic background, (ii) the 
attenuation of the neutron beam in the sample, and (iii) the 
magnetic form factor of Yb 3+ . The normalization is per for- 
mula unit (f.u.) of Yb 2 03. Data from runs with two incident 
neutron energies, 75 meV and 250 meV, are included. 



on flat samples whose thickness is small compared to 
their lateral dimensions (as is the case here) has been 
verified previously by Monte Carlo simulations incorpo- 
rating realistic model scattering cross-sections, 18 and also 
accords with experience from many other similar exper- 
iments. The background derived this way is included in 
Fig.n and is generally a little higher than the intensity 
from Y0.5LU0.5N12B2C. 

As mentioned earlier, the polycrystalline sample of 
YbN^b^C used in these experiments was contaminated 
with 9% (by mass) of Yb2C>3 impurity. To correct for the 
signal from this Yb2C>3 we measured the energy spectra 
of polycrystalline samples of pure Yb2C>3 and Y2O3, the 
latter playing the part of the non-magnetic reference sam- 
ple. We are not aware of any previous measurements of 
the magnetic excitations in Yb2 03, and so for reference 
we show in Fig.|3the energy spectrum of Yb2C>3 after cor- 
rection for the non-magnetic scattering by the method 
described above. Corrections have also been made for 
the attenuation of the neutron beam in the sample and 
for the free-ion magnetic form factor of Yb 3+ . Hence, the 
quantity plotted on Fig.^is 5(0, u), the zero-Q magnetic 
response function. Data from runs with two different en- 
ergies has been included. The strongest magnetic signal 
in the measured energy range is seen to be centred near 
70 meV. 

Figure^] shows the zero-Q magnetic response function 
of YbN^B^C derived from the runs with incident neu- 
tron energies 35 meV, 75 meV and 250 meV. To arrive at 
Fig.[3]we corrected the raw data for (i) the non-magnetic 



scattering shown in Fig.^ (ii) the Yb2C>3 impurity scat- 
tering, (iii) the attenuation of the neutron beam in the 
sample, and (iv) the magnetic form factor of Yb 3+ . The 
good agreement between the results from the three runs 
in the energy ranges where they overlap gives us confi- 
dence that the method used to estimate the non-magnetic 
signal is reliable. The Yb2 03 impurity correction has lit- 
tle effect (<10%) over most of the energy range, and 
is only significant around 70meV where it accounts for 
~40% of the raw magnetic signal in Fig.^ 

The dynamic magnetic response of YbN^B^C is seen 
from Fig. OH to be broad in energy extending beyond 
150 meV, but does show some structure. There are three 
peaks, the first centred close to zero energy, the second 
just below 20meV, and the third just above 40meV. 
These peaks can be seen more clearly in Fig.01 which 
highlights the low-energy data (up to 70meV) from the 
35meV and 75meV runs. 

To characterize the magnetic response we use Eq. (0J 
with a phenomenological spectral-weight function given 
by 



F(u) 



3 

EC-i 

i=l 



{hu + Eif 



(hco - Etf + r? 



( 5 ) 

i.e. the sum of three pairs of Lorentzian functions centred 
on ±2?i with energy widths (half-width at half maximum) 
IV The Cj coefficients satisfy Cj = 1 as required by the 
normalization condition, Eq. @. A Lorentzian function 
is a reasonable approximation when the peak broaden- 
ing is relatively small and arises from simple relaxational 
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FIG. 3: Magnetic response function of polycrystalline 
YbNi2B2C measured on the HET time-of-fiight spectrometer. 
The data are corrected for (i) the non-magnetic background, 
(ii) the YD2O3 impurity scattering, (iii) the attenuation of 
the neutron beam in the sample, and (iv) the magnetic form 
factor of Yb 3+ . The normalization is per formula unit (f.u.) 
of YbNi2B2C. Data from runs with three incident neutron en- 
ergies, 35meV, 75meV and 250 meV, are included. The solid 
line is calculated from Eqs. I^J an d with the parameters 
given in Table [I] 



processes. Later on we will compare the data with the 
response function calculated by an approximate method 
applicable to crystal held transitions in heavy fermion 
systems^ 

The solid lines drawn on Figs.0and|3|show that Eq. (J5J 
can provide a good description of the data. The pa- 
rameters used to fit the model response function to the 
data are given in Table [I] It is particularly satisfy- 
ing that the value Xau = 8.1 x 10~ 31 m 3 per ion (SI 
units), which corresponds to 0.039 emu/mol in Gaussian 
cgs units, matches closely with the powder susceptibility 
determined by magnetometry, which at 5 K is approx- 
imately 0.04emu/moliiS The justification for the first 
peak being quasielastic (i.e. centred on huj = 0) comes 
mainly from the low energy data to be presented shortly. 
The fit indicates that the 18meV and 43meV inelastic 
peaks are broader than the quasielastic peak. 

We now turn to the low energy part of the mag- 
netic response function. FigureEfa) displays the neu- 
tron inelastic scattering from the single-crystal samples 
of YbNi2B2C and LUM2B2C measured at a temperature 
of 1.9 K on the IN6 spectrometer. In these runs the crys- 
tals were oriented so that across the whole detector bank 
the scattering vector Q was approximately parallel to 
the a* reciprocal lattice vector. The two spectra are in- 
distinguishable in the negative energy region (neutron 
energy gain scattering), but differ appreciably at posi- 
tive energies (neutron energy loss) . The difference in the 



FIG. 4: Low energy (huj < 70meV) part of Fig.[3] showing the 
data from runs with incident energies 35meV and 75meV. 
The solid line is calculated from Eqs. @ and JKJ with the 
parameters given in Table and the broken lines depict the 
three component peaks that make up the overall line shape. 



TABLE I: Peak parameters of the Lorentzian spectral- weight 
function Eq. © used to describe the magnetic response func- 
tion of YbNi 2 B 2 C shown in Figs. [U and The value of the 
susceptibility determined from the fit is Xa« = 8.1 x 10 -31 m 3 
per ion in SI units, which corresponds to 0.039 emu/mol in 
Gaussian cgs units. 



Peak 


r " 


Ei (meV) 


Ti (meV) 


1 


0.91 ±0.03 


0.0 


2.1 ± 0.1 


2 


0.06 ±0.01 


18.2 ±0.7 


5.5 ± 1.1 


3 


0.03 ±0.005 


43.1 ± 0.9 


8.1 ± 1.5 



"The error given for c» is just the statistical error from the fit, and 
does not include the ~10% systematic uncertainty in the absolute 
calibration of the data. 

''The fitted width of peak 2 is from the quasielastic fit to the low 
energy data shown in Fig.lSl 



energy-loss scattering corresponds to the dynamic mag- 
netic response of the Yb ions in YDN12B2C. At negative 
energies the scattering from the two samples is essentially 
the same because the energy-gain magnetic scattering is 
strongly suppressed by the factor {l + n(u>)} in the cross- 
section, Eq. J2|. 

In order to isolate the magnetic scattering we must 
estimate the non-magnetic background. On IN6 the 
sources of background are different to those described 
earlier in connection with the HET measurements. Out- 
side of the energy range [—1, 1] meV the background is 
almost sample-independent at low temperatures, while 
between — 1 meV and 1 meV there is additional non- 
magnetic scattering from the elastic peak. This latter 
scattering scales with the amplitude of the elastic peak, 
which is sample-dependent. We confirmed these prop- 



6 



(a) Raw data 



3 b> 
E 




YbNi,B 9 C 



Energy (meV) 
(b) Background-corrected data 



3 b> 

O CD 

w E 

CO 

E 




Energy (meV) 

FIG. 5: (a) Neutron inelastic scattering from single crystal 
mosaic samples of YbNi2B2C and LUM2B2C measured on 
the IN6 spectrometer. The elastic peak is shown reduced by 
a factor 500 to indicate the energy resolution, (b) Dynamic 
magnetic response measured with the scattering vector ap- 
proximately parallel to the reciprocal lattice vectors a* and 
c*. The data have been corrected for (i) the non-magnetic 
background, (ii) the attenuation of the neutron beam in the 
sample, and (iii) the magnetic form factor of Yb 3+ . The solid 
lines show the model response function of Eqs. (|1J and 
Apart from an overall scale factor the parameters are the same 
as those listed in Table [I] 



erties of the background by comparing the LuNi 2 B 2 C 
spectrum with similar measurements from a vanadium 
sample and an empty aluminium sample holder. Our 
procedure to estimate the non-magnetic background was 
then as follows. First, we fitted a polynomial to the parts 
of the LuNi2B2C spectrum that lie outside [—1,1] meV. 
Next we subtracted this polynomial from the LuNi2B2C 
spectrum to separate the elastic peak, and scaled up this 
elastic peak to match that in the YbNi2B2C data. We 
added the scaled elastic peak back onto the polynomial, 
and substituted the result back into the LuNi2B2C spec- 
trum in place of the original [—1, 1] meV data. 
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FIG. 6: Dynamic magnetic response of YbNi2B2C in the 
low energy range up to 6meV. Data from three different 
measurements are shown, each separately corrected for the 
non-magnetic background. The lowest-energy data are the 
Q ~|| a* points shown on Fig.^b), and the other runs are 
scaled to match in the energy range where they overlap. The 
solid line is the model response function of Eqs. gj and JSJ 
calculated with the parameters listed in Table apart from 
ci which was adjusted to fit the present data. 



In Fig.[S{b) we show the magnetic response for two 
crystal orientations after subtraction of the modified 
LuNi 2 B 2 C spectrum from the YbNi 2 B 2 C spectrum. Cor- 
rections have also been made for the attenuation of the 
neutron beam in the sample, and for the magnetic form 
factor of Yb 3+ . The magnetic response when Q is ap- 
proximately parallel to a* is seen to be a factor two larger 
than when Q is approximately parallel to c*. Measure- 
ments on the triple-axis spectrometer IN14, which in- 
volves averaging over a much smaller range of Q than on 
IN6, were also consistent with a factor two difference in 
relative intensities for these orientations. 

The energy range covered by the data Fig.[5[b) is too 
small to allow us to examine the line shape of the low 
energy magnetic response in any detail, and so in Fig-IHl 
we have plotted the Q ~|| a* points from Fig.^b) to- 
gether with data from two other runs appended so as 
to extend the energy range up to 6 meV. The additional 
measurements were made on IN6 with an incident neu- 
tron energy of 4.9 me V and on IN 14 with a fixed final 
energy of 4.7 meV. The same single-crystal samples were 
used for all the measurements, and after correction for 
the non-magnetic background the new data points were 
scaled so that all three runs matched up in the energy 
ranges over which they overlap. 

We used the data shown in Fig.|f)|to refine the model 
response function defined by Eqs. and JSJ. In this low 
energy range the 18meV and 43meV peaks have negli- 
gible weight, as can be seen from Fig.0J and the model 
function is dominated by the first peak. Allowing the 
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parameters of peak 1 to vary we achieved the best agree- 
ment with a quasielastic peak (i.e. ~E\ — 0) with width 
Ti = 2.1 ± 0.1 meV. Attempts to force the peak to be 
inelastic (E\ > 0) resulted in progressively poorer fits as 
the centre of the Lorentzian was shifted to higher ener- 
gies. These values of E\ and Ti were then used as fixed 
parameters in a further fit to the polycrystalline data 
(Figs -El and 0J) which finally established the other model 
parameters listed in Tabled We note that a powder aver- 
age of the single-crystal data shown in Fig.|Sf b) agrees to 
within 25% of the the model response function displayed 
on Figs. [31 and 0] which is an acceptable margin of error 
given the uncertainties in latter fit and in the absolute 
calibration of the data. This agreement gives us further 
confidence in the Ci and Xav amplitude parameters. 

We mention finally that on the IN14 triple-axis spec- 
trometer we examined the low-energy scattering in a se- 
ries of energy scans at fixed Q, with Q chosen at differ- 
ent points along the main symmetry directions (1,0,0), 
(1, 1, 0) and (0, 0, 1), as well as at some off-symmetry po- 
sitions. As already mentioned, the intensity of the scat- 
tering was found to vary with the direction of Q (smallest 
when Q is parallel to the c* direction and constant within 
the a*b* plane), but to within the statistical precision of 
the data we did not observe any change in the line shape 
with Q. The measurements allow us to place an upper 
limit of 0.5 meV on the extent of any dispersion in the 
low energy magnetic response. 



IV. ANALYSIS AND DISCUSSION 

These measurements have revealed that the magnetic 
excitation spectrum of YbNi2B2C contains three peaks, 
centred approximately on OmeV, 18meV and 43meV, 
with widths of several meV. The broadening is such that 
there is significant overlap of the peaks, but not enough 
to prevent the peaks from being resolved. 

What does this tell us about the Yb electronic state in 
YbNi2B2C? In general, the properties of Yb intermetal- 
lic compounds are influenced by a partial hybridization 
of the Yb 4/ electrons with conduction or valence elec- 
tron states of the hosti^S When the hybridization is very 
weak the Yb ions are trivalent with localized 4/ electrons 
carrying a magnetic moment. The dominant external in- 
teraction is then with the crystalline electric field (CEF), 
which splits the 2 F 7 / 2 ground state term of Yb 3+ into a 
manifold of levels. Transitions between these CEF levels 
give rise to sharp peaks in the magnetic excitation spec- 
trum. The opposite extreme is when the hybridization is 
much stronger than the CEF. This results in valence fluc- 
tuations of the Yb ions, and the corresponding dynamic 
magnetic response is very broad in energy. 

The case of YbNi 2 B 2 C is seemingly part-way between 
these extremes. However, the fact that we observe in- 
elastic peaks rather than a single broad response sug- 
gests that a description in terms of localized 4/ electrons 
in a CEF is the better starting point. This scenario is 



consistent with the X-ray absorption spectrum measured 
at the Lni edge, which indicated a stable 3+ ionization 
state for the Yb ionsi 

Given a stable Yb valence, the broadening of the CEF 
transitions is expected to arise from the same strong spin 
fluctuations that dominate the electronic heat capacity at 
low temperatures. To see if this is the case we consider 
the relation 



obtained by Edwards and Lonzarich 21 for the low- 
temperature electronic heat capacity 7T of a system of 
fluctuating spins described by a set of over-damped har- 
monic oscillators with a Q-independent relaxation width 
r. Here, v is the number of fluctuating degrees of free- 
dom. Eq. © has been shown to hold for a range of 
materials in which strong spin fluctuations are known to 
be important^ For an effective spin-^ system, such as 
we have here, the spin state is defined by one component 
and so v = 1. If we take T = 2.1 meV, the quasielas- 
tic width, then Eq. © predicts 7 = 0.36 J mol -1 K~ 2 , 
which is roughly 30% smaller than found experimentally 
for YbNi2B2C. This reasonable level of agreement sup- 
ports the hypothesis that spin fluctuations account for 
most of the low temperature electronic heat capacity in 
YbNi 2 B 2 C. Exchange interactions between the Yb ions 
can be ruled out as a significant contributor to the CEF 
broadening because the magnetic excitations do not ex- 
hibit any observable dispersion. Furthermore, if Yb-Yb 
interactions were important then dilution of the Yb ions 
by non-magnetic Lu should make the magnetic excita- 
tions sharper and therefore easier to observe, but neu- 
tron measurements on Ybo.iLuo.gNi2B2C have failed to 
detect any such sharpening^ 

We will now examine the CEF in YbNi2B2C more 
quantitatively. CEF energy spectra for i?Ni2B2C (R = 
Dy, Ho, Er and Tm) measured by neutron spectroscopy 
have been presented and analyzed in detail by Gasser 
et alm^SA For the 4/mmm point symmetry of the 

R site the angular dependence of the CEF is described 
by five parameters (essentially coefficients of spherical 
harmonics). Gasser et al. derived values for the CEF 
parameters for R = Dy, Ho, Er and Tm by comparing 
quantities calculated from the CEF model with the avail- 
able spectroscopic and magnetic data. After correction 
for the expected dependence on the ion size there is little 
variation in the CEF parameters for each ion. Hence, 
Gasser et al. extrapolated the CEF parameters to other 
R ions and predicted the CEF splittings. 

For the case R = Yb, the extrapolations predict that 
the 2 F 7 / 2 term splits into four Kramers' doublets. The 
predicted energies and wavefunctions of the doublets are 
given in Table [H] The neutron scattering cross-sections 
for transitions (i) within the ground state, and (ii) from 
the ground state to each of the three excited doublets 
are calculated to be (i) 497mbsr~ 1 [f.u.] -1 , and (ii) 260, 
228 and 9mbsr _1 [f.u.]" 1 . From Fig. 0]it is evident on 
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TABLE II: Energy levels, symmetries and wavefunctions of 
the CEF doublets in YbNi 2 B 2 C calculated from the CEF ex- 
trapolated from TmNi2B2Ci^ The CEF parameters (Stevens' 
operator notation) are B° = -3.95 x 10" \ = -3.46 xlO" 3 , 
B\ = 1.01 xlO" 1 , Be = -1.77 xl0~ 4 , fl| = 2.09 x 10" 3 meV. 



Energy 


Symmetry 


Wavefunction 6 


(meV) 


label" 




0.0 




0.902| ± 7/2) - 0.431| =F 1/2) 


7.9 


r-7 


0.717| + 3/2) - 0.697| ±5/2) 


24.3 


r 6 " 


0.431 ± 7/2) + 0.902 + 1/2) 


26.0 


17 


0.697 +3/2) +0.717 ±5/2) 



"Irreducible rep rese ntations of the A/mmm (D^h) double group, 
as defined in Ref . l25l 

b \mj) is shorthand for |J, mj), where J = 7/2 is the combined 
spin and orbital angular momentum quantum number of the 2 Ff/2 
term, and raj is the magnetic quantum number. We assume neg- 
ligible admixture of states from the 2 F 5/ /2 term. In the iJNi2B2C 
structure the CEF quantization axis is parallel to the crystallo- 
graphic c axis. 

the one hand that the observed CEF splitting is nearly 
twice that predicted by the extrapolated CEF model, but 
on the other that the relative spacings of energy lev- 
els and their transition intensities match quite closely 
with the predictions. We see also that the absence of a 
third excited CEF peak in the experimental spectrum is 
explained by its small spectral weight and by its close 
proximity to the second excited CEF peak. 

Encouraged, we now consider the anisotropy in the 
magnetic scattering (Fig. 0J>) and bulk magnetic sus- 
ceptibility. The neutron scattering cross-section for a 
transition between two CEF levels % and j depends on 

Ja\j)\ 2 , the squared matrix elements of the angular 
momentum operator J a (a — x,y,z)*& The (1 — Q 2 ) 
factor in Eq. [21 means that when Q || c* the intensity is 
proportional to J x \j)\ 2 ± \{i\Jy\j)\ 2 1 and when Q || a* 
the intensity is proportional to J y \j)\ 2 + Kil^lj)! 2 - 

In tetragonal symmetry |(«|</z|j)| 2 = \(i\J y \j)\ 2 , and 
assuming the low-energy scattering comes mainly from 
scattering within the ground-state doublet we deduce 
from the relative intensities in Fig. |3Jd that |(0|J Z |0)| 2 » 
3| (0| J x \0) | 2 . However, from the wavefunctions of the ex- 
trapolated CEF model (Table OH) we find |(0|J Z |0)| 2 « 
50 1 (0| J x \0) | 2 - In other words, the sense of the anisotropy 
predicted by the extrapolation is correct, but the ob- 
served degree of anisotropy is much smaller than pre- 
dicted by this CEF-only model. 

The same conclusion can be drawn from the bulk sus- 
ceptibility. Figure compares the susceptibility calcu- 
lated from the extrapolated CEF model with the sin- 
gle crystal data reported in Ref. 0. Again, the sense of 
anisotropy is in accord, but the model predicts a much 
greater degree of anisotropy than observed experimen- 
tally. A further observation is that the measured powder- 
averaged susceptibility x&v is less than the calculated x&v , 
most of the difference being in the H \\ c component. At- 




Temperature (K) 

FIG. 7: A comparison of the measured magnetic susceptibility 
of YDN12B2C (open circles) with the susceptibility calculated 
from the scaled CEF model (lines). The experimental data 
has been taken from Avila at aL- The inset shows the inverse 
of the powder-averaged susceptibility. 

tempts to refine the CEF model indicated that a reduc- 
tion in magnetic anisotropy could be achieved only at the 
expense of a poorer agreement with the scattering inten- 
sities and without significant change in x&v- Hence, we 
were not able to find an acceptable alternative model that 
was substantially different from the extrapolated CEF 
model given in Table ITT1 

The discrepancy in x&v increases dramatically at low 
temperatures, but when compared on a plot of x&v vs 
T, as done in the inset to Fig. [7| it becomes apparent 
that over much of the temperature range the measured 
and predicted curves have the same slope. This suggests 
the influence of an effective antiferromagnetic exchange 
interaction on the Yb ions, a likely source of which is 
the on-site Kondo interaction (we have already ruled out 
the existence of measurable inter-site antiferromagnetic 
interactions). In the following analysis, therefore, we con- 
sider the addition of a Kondo interaction to the CEF po- 
tential established above, and use approximate methods 
to see whether this model leads to a better description of 
the experimental results. 

We consider first the influence of a Kondo interaction 
on the dynamic magnetic response. Previously, a sat- 
isfactory description of the measured neutron magnetic 
scattering spectra of several Yb Kondo lattice compounds 
has been obtained from a simple approximate solution 
to the Anderson impurity Hamiltonianii242& Within this 
scheme (the ZZF approximation^) one can calculate the 
single-ion dynamic magnetic response using as input the 
CEF parameters, the Kondo temperature 7k, and the 
number of / electron holes rif . To apply this method to 
YbNi2B2C we scaled the CEF parameters extrapolated 
from TmNi2B2C by a factor of 1.8 to match the energy 
scale of the observed CEF splitting, and set rif = 0.95, 
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FIG. 8: Dynamic magnetic response of YbNi2B2C calculated 
by the method of Ref. compared with the experimental 
data from Fig. 0] The model contains the CEF parameters 
extrapolated from TmNi2B2C scaled by a factor of 1.8, and 
a Kondo temperature Ik = 25 K. 

close to unity to reflect that the Yb valence is close to 3+. 
When 0.9 < n./ < 1 the results are not sensitive to the 
precise value chosen for rif. We fixed Xk = 25 K from 
the relation fceTk = T = 2.1meV from the quasielas- 
tic width. 27 Fig. [SJshows the response calculated by the 
ZZF approximation together with the HET data from 
Fig. Given that no attempt has been made to re- 
fine the CEF model the level of agreement as far as the 
peak positions and relative intensities are concerned is 
surprisingly good. This indicates that the scaled CEF 
model is a good starting point for a description of the 
CEF in YbNi 2 B 2 C. What is also clear from Figgis that 
the ZZF approximation appreciably underestimates the 
widths of the 18meV and 43meV peaks. A complete 
understanding of the lineshape is evidently lacking 2 ^ 

We next turn to the powder-averaged magnetic sus- 
ceptibility, which is connected to the dynamic magnetic 
response by Eq. Integrating the ZZF dynamic re- 
sponse given in Fig.|S|we obtain Xav = 0.036 emu/mol (or 
1 /Xav = 28 mol/emu) at X = 5 K. This result is very close 
to the experimental value, as may be seen from the inset 
to Fig. [3 At higher temperatures, however, we find that 
the ZZF model systematically exceeds the experimental 
susceptibility as temperature increases. In order to ob- 
tain agreement it is necessary to increase Xk in the ZZF 
model systematically with temperature. We find that the 
Tk needed in the ZZF model to achieve agreement with 
experiment increases smoothly with temperature, even- 
tually saturating at Xk = 130 K at high temperatures. 

Such a renormalization of the Kondo interaction with 
temperature is well known in systems where the CEF 
splitting is larger than the low temperature Xk, and has 
already been discussed for YbNi 2 B 2 C by Rams et alm^ in 
connection with measurements of the electric field gradi- 
ent. The effective Kondo interaction renormalizes down 



with decreasing temperature as a result of the reduc- 
tion in the effective degeneracy of the Yb 3+ ion due to 
the CEF splitting. According to theory^ the effective 
Kondo temperature at high temperature X^ is related to 
that at low temperature X^ by 

X£ « (AxA^Xk) 1 / 4 , (7) 

where the Ai are the energies (in Kelvin) of the excited 
CEF levels relative to the ground state. Eq. Q) is valid 
if X^ < A„ a condition that is satisfied here. If we set 
X K = 25 K and use the observed CEF levels for Aj then 
Eq. Q predicts X^ 1 ~ 190 K, which is comparable to the 
value Xk = 130 K required in the ZZF model to match 
the observed high temperature susceptibility. 

Finally, we comment on the neutron scattering mea- 
surements reported earlier by Sierks et alM> In Ref. ^| 
the authors investigated both the low and high energy 
part of the excitation spectrum using a polycrystalline 
sample of YbNi 2 B 2 C. Without a non-magnetic reference 
spectrum they were not able to find the broad inelastic 
magnetic peaks centred on 18meV and 43meV we have 
identified here. They did, however, observe two inelastic 
magnetic features at low energies. One of these is a broad 
shoulder on the side of the elastic peak extending up to 
approximately 1 meV, which developed into a quasielas- 
tic peak with increasing temperature, and the other a 
broader peak centred near 3.5 meV. We believe the lat- 
ter feature is what we have identified here as quasielastic 
scattering within the ground state doublet, but concern- 
ing the scattering below 1 meV there is no sign of this 
in our data, as is clear from Fig. EJd. The possibility 
that this scattering is absent from our data because it 
is highly localized in reciprocal space can be discounted 
because the large detector bank of the IN6 spectrometer 
used to collect the single-crystal data meant that the two 
scans shown in Fig. 03 include data from all parts of the 
Brillouin zone. 

The most likely explanation, therefore, is that the 
magnetic scattering observed by Sierks et al. below 
1 meV originates from an impurity in their sample, 
probably Yb 2 03. To support this argument we show 
in Fig. the background-corrected spectrum of our 
own polycrystalline sample of YbNi 2 B 2 C and in Fig. 

the background-corrected spectrum of polycrystalline 
Yb 2 03, both spectra measured on the IN5 spectrometer 
under the same experimental conditions to those used by 
Sierks et al. The peak centred at 0.5 meV present in both 
spectra is clearly a signature of Yb 2 03j2i and its inten- 
sity ratio in Figs.|5Ja, and b (when scaled by the respective 
formula masses) is consistent with the 9% mass fraction 
of Yb 2 03 found in the phase analysis of our YbNi 2 B 2 C 
sample. Assuming the same impurity was responsible for 
the 0.5 meV signal in Sierks et al's data we estimate that 
their YbNi 2 B 2 C sample contained ~4% Yb 2 3 . 
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FIG. 9: Low-energy magnetic excitations measured by neu- 
tron inelastic scattering from polycrystalline samples of (a) 
YbNi2B2C and (b) Yb203. The data have been corrected for 
(i) the non-magnetic background, (ii) the attenuation of the 
neutron beam, and (iii) the magnetic form factor of Yb 3+ . 
The data show that the peak centred near 0.5 meV in the 
YbNi2B2C data originates from the ~9% Yb203 impurity 
present in the sample. 



V. CONCLUSIONS 

The results and analysis presented here reveal that the 
magnetic properties of YDM2B2C are principally deter- 
mined by CEF and Kondo interactions acting on the 
Yb 3+ ions, with the CEF about an order of magnitude 
larger than the Kondo energy scale (Ik ~ 10 K from the 
single-impurity expressions for the heat capacity and sus- 
ceptibility, and Tk ~ 25 K from the quasielastic width) 2 ' . 
Any magnetic coupling between Yb ions is negligible. 
The CEF potential is roughly a factor 2 greater than 
found in other .RN12B2C compounds, presumably ampli- 
fied by hybridization of the Yb 4/ electrons with neigh- 
bouring atomic orbitals. Being an electrostatic interac- 
tion, hybridization is expected to modify the CEF to 



some degree; here, the main effect seems to be a uniform 
enhancement of the CEF without significant change in 
its anisotropy. 

The ZZF approximation scheme for the Anderson im- 
purity Hamiltonian succeeds in connecting the low tem- 
perature magnetic susceptibility and dynamic magnetic 
response, and thus provides a microscopic basis on which 
to interpret the properties of YbNi2B2C. The effective 
Kondo scale increases with temperature due to the ther- 
mal population of excited CEF states and associated in- 
crease in orbital degeneracy and effective moment. This 
mechanism explains why the Weiss temperature 9 p de- 
rived from the high temperature susceptibility is much 
larger than the value of the Kondo temperature deduced 
from various low-temperature properties. In effect, 9 P 
is a measure of the high temperature Kondo scale. It 
would be interesting to perform inelastic neutron scat- 
tering measurements at higher temperatures to test the 
ZZF prediction of a marked broadening of the dynamic 
magnetic response with increasing Tk. 

Our work has shown that in most respects YbNi2B2C 
can be understood in terms of conventional ideas for 
Kondo-lattice intermetallic compounds. However, in one 
respect, not mentioned so far, there are indications that 
less conventional physics may be needed. This concerns 
the low temperature thermodynamic, transport and mag- 
netic data which, if one examines the literature closely, 
are seen to exhibit significant deviations from Fermi liq- 
uid behaviour. Specifically, the susceptibility does not 
saturate below Tkr*^ the sommerfeld coefficient 7 con- 
tinues to increase as T tends to zero^IL and the low tem- 
perature resistivity could be consistent with a T n depen- 
dence with n < 2& 

These temperature dependences in fact bear a strong 
rescmblence to non-Fermi-liquid effects exhibited by 
YbRh2Si2, a material recently proposed to be the first 
example of a clean Yb compound showing quantum- 
critical behaviour above a low-lying antiferromagnetic 
transition^ We conjecture, therefore, that at zero tem- 
perature YbNi2B2C may be close to a quantum-critical 
point on the non-magnetic side. The ground state would 
then contain strong quantum fluctuations of the spins, 
and deviations from the Kondo model should be ex- 
pected. Indeed, the fact that Eq. ©, which is based on 
a phenomenological model for the spin fluctuation spec- 
trum, succeeds in connecting the experimental values of 
7 and r (to within 30%), whereas the single- impurity 
Kondo model does not (the Tk values deduced from 7 
and r differ by a factor of 2.5) may be evidence for a 
breakdown of the conventional heavy fermion scenario. 
Unfortunately, the present experiments are not sensi- 
tive enough to shed any further light on this possibility, 
but the indications are clear enough to suggest that fu- 
ture investigations into the low temperature properties 
of YbNi2B2C may prove to be fruitful. 
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